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We study the dynamics in the neighborhood of fixed points in a 4D symplectic map by means of 
the color and rotation method. We compare the results with the corresponding cases encountered 
in galactic type potentials and we find that they are in good agreement. The fact that the 4D 
phase space close to fixed points is similar to the 4D representations of the surfaces of section 
close to periodic orbits, indicates an archetypical 4D pattern for each kind of (in)stability, not 
only in 3D autonomous Hamiltonian systems with galactic type potentials but for a larger class 
1^ ■ of dynamical systems. This pattern is successfully visualized with the method we use in the 

O ■ paper. 
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1. Introduction 

Patsis and Zachilas [1994] proposed a method to visualize 4D spaces of section. It is based on rotation of 
the 3D projections of the figures in order to understand the geometry of the projections and on color for 
understanding the distribution of the consequents in the 4th dimension. This combined information from 
the rotated, colored, 3D projections allows us to associate specific structures with the dynamical behavior 
in the neighborhood of periodic orbits. 

Recently the method has been applied successfully in a series of papers that studied the dynamics in 
the neighborhood of periodic orbits in a 3D galactic potential [Katsanikas and Patsis 2011, Katsanikas 
et al 2011a, Katsanikas et al 2011b]. The results of these papers, together with those of the original one 
by Patsis and Zachilas [1994], as well as further results from work in progress, find a consistency between 
the shapes of the encountered structures in the 4D spaces of section and the kind of (in)stability of the 
periodic orbit. Despite the fact that until now no exceptional behavior has been found, the results cannot be 
characterized as generic. The up to now studied systems are 3D autonomous Hamiltonians with potentials 
suitable to describe stellar motion in 3D rotating galactic disks. They can be used to investigate observed 
morphological features out of the equatorial planes of disk galaxies [Patsis et al 2002] . 
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The motivation for investigating the phase space structure of a 4D symplectic map is to check whether 
or not the behavior encountered in the Hamiltonian systems is model dependent. So, we apply the same 
technique in the case of a 4D symplectic map and we study the structure of the phase space in the case of 
stability and in cases of instabilities studied in the Katsanikas' papers. Our goal is to compare the dynamics 
of the 4D map with those found in the Hamiltonian system, testing in this way the ubiquity of the results 
of the latter studies. 

The paper is structured as follows: In Sec. [2] we briefly describe the method, in Sec. [3] we present the 
map we use in our study, in Sec. H] we describe our results and finally we enumerate our conclusions in 
Sec. El 



2. A brief description of the method 

We consider the map T : R 4 — >• R 4 and follow the evolution of the array (x 4 , x 2 , X3, x 4 ). A set of three 
coordinates, e.g. (x±, X2, X3), are used for the 3D projection, while the fourth coordinate (e.g. x 4 ) determines 
the color of the consequents. There is a normalization of the color values in the [min(x 4 ), max(x 4 )] interval, 
which is mapped to [0,1]. In order to understand the shape of each 3D projection we rotate the figures 
on the screen of our computer. For presentations on paper we chose a set of projection angles that help 
the reader understand the shape of the object we describe in the text. We use in our applications the 
"Mathematica" package. Following the intrinsic "Mathematica" subroutines our viewpoint is given in 
spherical coordinates. The unit for the distance d of the consequents of the surface of section from the 
observer is given by "Mathematica" in a special scaled coordinate system, in which the longest side of the 
bounding box has length 1. For all figures we use d = 1. The method associates the smooth distribution or 
the mixing of colors, with specific types of dynamical behavior in the 4th dimension [Patsis and Zachilas 
1994, Katsanikas and Patsis 2011, Katsanikas et al 2011a, Katsanikas et al 2011b]. For a more detailed 
description see [Katsanikas and Patsis 2011]. 



3. The map 

The map we have chosen, T, belongs to a family of nonlinear symplectic 4D mappings in R 4 that is a 
generalization of the standard map. The definition in the general form is: 



/ x\\ / x\ + K\ sin(xi + x 2 ) + L\ sin(xi + x 2 + x 3 + x 4 ) \ 
x 2 _ X\+ X2 

x 3 + K 2 sin(x 3 + x 4 ) + L 2 sin(xi + x 2 + x 3 + x 4 ' 
\x 3 + x 4 



£3 

\x 4 y 



(mod 2ir) 



Several cases of the maps of this family have been used in the past to study the dynamics in the 
neighborhood of fixed points. Already Froeschle [1972] used (1) with L\ = L 2 and tried to visualize the 
4D surfaces of section. This was an additional motivation for choosing the particular map to apply our 
visualization technique. Furthermore studies accomplished by Pfenniger [1985], Olle and Pfenniger [1999], 
Jorba and Olle [2004], guarantee that Hamiltonian Hopf bifurcations are happening in this system. This 
allows us to compare also the behavior of the galactic type Hamiltonian with that of the map at transitions 
from stability to complex instability. 

Following Pfenniger [1985] and Jorba and Olle [2004], we examine the case with L\ + L 2 = having as 
conjugated variables (x 4 ,x 2 ) and (x3,x 4 ). We also restrict the parameter space by taking K 2 = 0. Then, 
we rename L\ = —L 2 as L, and K\ as K. It is evident from the stability diagrams in Pfenniger [1985], that 
there is always a critical value L c = — if/4, for —8 < K < 0, for which we have a transition from stability 
to complex instability. In order to include the case of such a transition in our calculations, we take here 
K = -2. 

We performed our calculations in double, as well as in quadruple precision and we found the same 
results, except in a case of simple instability (see Sect. |4] below), where the results are only qualitatively 
similar. In that case we present the results found with quadruple precision. 
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Fig. 1. The orbital behavior in the neighborhood of a stable fixed point in the L = 0.3 case, (a) The characteristic toroidal 
figure associated with stability. The arrows indicate the regions where the smooth color variation changes from the external 
to the internal side of the torus. The point of view is (0,<j>) = (10°, 270°) (b) An enlargement of (a) at one of the transition 
regions of the color from the external ( u Ext") to the internal ("Int") side of the torus (cf. with figure 12 in Katsanikas and 
Patsis [2011]). We use less consequents in order to show how the structure of the colors varies. The point of view is now 
(8, (f>) = (18°, 40°) to show better how the color change side. 



4. 4D phase space 

We studied the dynamics close to the fixed point (0,0,0,0) by varying L. In the neighborhood of stable 
fixed points < L < 0.5 we find in the 4D phase space invariant tori. Such invariant tori have been 
presented in projections already by Froeschle [1972]. However, with the method of color and rotation we 
could verify that these tori are 4D objects. Their 3D projections are of the kind called "rotational tori" 
by Vrahatis et al. [1997]. Their 4D representation, in all cases we examined, is in perfect agreement with 
the rotational invariant tori with smooth color variation on their surfaces found by Katsanikas and Patsis 
[2011] (see their figure 13). In Fig. [1] we give a typical example for L = 0.3. The initial conditions are 
(xi, X2, xs, X4) = (0.2,0,0,0). We use the (xi,X2,xs) 3D projection for the spatial representation of the 
consequents and X4 to give the color that represents the fourth dimension. In this and in all subsequent 
figures the color bar on the right hand side gives the range of the color variation that corresponds to the 
normalized values of the fourth coordinate. Even the change of the smooth color variation from the external 
to the internal side of the torus described by Katsanikas and Patsis [2011] in their study of the autonomous 
Hamiltonian system is observed in the case of the rotational torus around the stable fixed point of the 
symplectic map. The detailed description of the color variation on the surface of the rotational torus should 
be sought in Katsanikas and Patsis [2011]. These structures characterize the dynamics of ordered behavior 
in the case of our map. Calculating Lyapunov Characteristic Numbers for all these orbits on the tori we 
find that they tend very fast to zero. 

Simple instability is encountered for —4 < L < as the |L| value increases. In the neighborhood of 
simple unstable (0,0,0,0) fixed points (two eigenvalues on the unit circle and two on the real axis) we 
find consequents defining a ribbon of the form of a double loop. This characteristic double loop structure 
is encountered only close to a critical value of the parameter (the value of L) for which the stability of 
the fixed point changes and it becomes simple unstable. A typical example for L = —0.03 and initial 
conditions {x\, X2, X3, X4) = (0.001,0,0,-0.001) is presented in Fig. [2l It gives the first 10 5 consequents. 
Such structures are associated with simple instability in 3D galactic potentials by Patsis and Zachilas 
[1994] and are extensively investigated by Katsanikas et al. [2011c]. The current situation resembles the 
dynamical behavior encountered in the neighborhood of simple unstable periodic orbits of the z-axis family 
in 3D rotating Hamiltonian systems [Katsanikas et al. 2011c], when we have a surface of section defined by 
z = 0. In the intersection of two branches of the double loop structure we have the same color, indicating 
that it belongs to the cases where the double loops correspond to a real 4D "8" -like structure. 
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Fig. 2. The orbital behavior in the neighborhood of a simple unstable fixed point in the L = —0.03 case. We observe the 
characteristic double loop figure associated with simple instability. The point of view is (6,(f>) = (120°, 9°). 
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Fig. 3. (a) The (2:3, £4) projection of the phase space close to a simple unstable fixed point for L — —0.1. The consequents 
depart from the double loop structure after about 14 x 10 iterations, (b) The LCN(n) evolution. We observe that LCN levels 
off at a positive value about 0.355. 

After a large number of iterations the consequents depart from the surface of the ribbon-like object 
and build a cloud of points filling finally all available phase space. As \L\ increases the consequents leave 
the "8" -like structure faster. In the example of Fig. [3j L = —0.1 and the initial conditions are as in the 
previous case. We observe that the consequents have left the "8" -like structure after 14 x 10 4 iterations and 
start visiting larger volumes of the phase space (Fig. [3h,). This is a case where we encountered quantitative 
differences when we use quadruple precision in our calculations. However, also in this case the results are 
qualitative the same. This behavior is depicted in the evolution of the "finite time" Lyapunov Characteristic 
Number LCN(n) = (1/n) In |£(n)/£(0)|, where £(0) and £(n) are the distances of two nearby points at 
n = and after n iterations [see e.g. Skokos 2010]. The LCN(n) for the example of Fig. [3k, is given in 
Fig. [3b. The evolution of the LCN(n) index has been followed for 10 6 iterations. The index has an overall 
decreasing part that ends after about 14 x 10 4 iterations (when the consequents start departing from the 
4D "8" -structure). Then it increases and finally it levels off at a value about 0.355. 

The fixed point (0,0,0,0) becomes double unstable for L < —4. In all cases we have examined, close to 
double unstable fixed points we observe clouds of scattered consequents with mixing of colors, even close to 
the transitions from simple to double instability. In the example of Fig. 0^ we give an orbit for L = —4.01 
with initial conditions (x\,X2, £3, X4) = (0.1,0,0,0). In this case LCN(n) levels off almost immediately at 
a value about 0.85 (Fig. HJd). This is the usual behavior also in most cases studied by Katsanikas et al 
[2011c] close to double unstable periodic orbits. 

For K = —2, the eigenvalues of the characteristic polynomial of the monodromy matrix become 
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Fig. 5. The orbital behavior in the neighborhood of a complex unstable fixed point in the L = 0.51 case. We observe a disky 
structure with smooth color variation in the 4D phase space of our symplectic map. The point of view is (6, <j>) = (3°, 216°). 

complex, off the unit circle, for L > 0.5. The critical value at which happens the transition from stability 
to complex instability is L = 0.5, i.e. for this value we have a Hamiltonian Hopf bifurcation. Studying 
the dynamics close to the (0,0,0,0) fixed point, for L = 0.51, just beyond the critical value, we find the 
disky structure (confined torus) that characterizes the dynamics close to complex unstable periodic orbits 
or fixed points. This structure has been found in the 3D projections of the 4D phase space by Pfenniger 
[1985], Jorba and Olle [2004] and Olle et al. [2004]. It has been shown by Katsanikas et al. [2011a] that this 
structure is a 4D object with a smooth color variation as we move from one side of the "disk" to the other. 
For the case of the map (1) we can see it in Fig. [5] for 2 x 10 5 iterations. These "disks" have an internal 
spiral structure [Contopoulos et al. 1994, Papadaki et al. 1995, Contopoulos 2002], the 4D structure of 
which has been shown by Katsanikas et al [2011a]. A similar internal spiral structure as the one described 
by Katsanikas et al [2011a] in the 4D confined tori (see their figures 2,3,4 and 5) has been encountered 
in all cases of map (1) we have studied in this paper. In the present case the spiral has six arms. The 
behavior of the LCN indices in the two cases are also similar. Fig. [6] depicts the variation of LCN(n) up to 
n = 10 5 and is similar to the LCN(t) variation of the corresponding case of the autonomous Hamiltonian 
system (cf. figure 6 in Katsanikas et al 2011a). In both case an initial decrease of the indices is followed by 
a leveling off. 

Also the behavior in the neighborhood of a complex unstable fixed point away from the critical L c value 
resembles that in the neighborhood of a complex unstable periodic orbit away from the Jacobi constant 
value for which we have the transition from stability to complex instability. We have clouds of points with 
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Fig. 6. The LCN(n) evolution in the case of Fig.[S] We observe that LCN(n) levels off at a positive value about 0.009, after 
an initial decreasing part. 

mixed colors. In the latter case the LCN(n) index reaches a positive value, larger than that of the confined 
torus case. The whole variation of the index remains in the cases we examined larger than the values we 
obtained for the disky structure for the same number of iterations. 

5. Conclusions 

In this paper we compared the dynamical behavior in the neighborhood of fixed points in a 4D symplectic 
map and we compared the results with those found in the neighborhood of periodic orbits in a 3D au- 
tonomous Hamiltonian system with a potential of galactic type. In our study we used the method of color 
and rotation. We found that the structures encountered in the spaces of section in the Hamiltonian system, 
which characterize the kind of (in) stability of the periodic orbits are also found characterizing the dynamics 
in the neighborhood of fixed points in the 4D symplectic map. These structures are: (1) The tori with the 
smooth color variation for stability, (2) the double loop ribbons with smooth color variation for the simple 
unstable case close to transitions from stability to simple instability, (3) the clouds of scattered points 
in the neighborhood of double unstable periodic orbits, and (4) the disky structures with smooth color 
variation close to a complex unstable fixed point near the transition from stability to complex instability. 
These results indicate that the 4D visualization of the surfaces of section in 3D autonomous Hamiltonian 
systems and of the phase space in 4D symplectic maps associates the dynamics in the neighborhood of a 
fixed point, or a periodic orbit, with specific structures. We have not found any exception even in the case of 
the symplectic map we used in the present study. Thus, these structures seem to reflect a generic behavior 
in a larger class of dynamical systems. Our study also indicates that the color and rotation method allows 
us to know the kind of instability of the fixed point by a simple inspection of the 4D phase space and gives 
in a direct way an insight of the dynamics in such regions, not only in galactic type potentials, but in a 
broader spectrum of dynamical systems. 
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